Transverse magnetic heat transport on the topological surface 
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We investigate Nernst-Ettingshausen and thermal Hall effects on the surface of a topological 
insulator on which a ferromagnet is attached. We present general expressions of the Peltier and the 
thermal Hall conductivities, which are reduced to simple forms at low temperatures. It is shown that 
the Peltier and the thermal Hall conductivities show non-monotonous dependence on temperature. 
At low temperature, they have linear dependence on temperature. From the behavior of the Peltier 
conductivity at low temperature, one can estimate the magnitude of the gap induced by time-reversal 
symmetry breaking. Moreover, we find that the Peltier conductivity can be used to map the Berry 
phase structure. These results open up possibility to control transverse heat transport magnetically. 

PACS numbers: 73.43.Nq, 72.25.Dc, 85.75.-d 



Topological insulator provides us with a new state of 
matter topologically distinct from the conventional band 
insulator^. In particular, edge channels or surface states 
are described by Dirac fermions and protected by the 
band gap in bulk states. As a result, backward scatter- 
ing is forbidden by time-reversal symmetry. From the 
viewpoint of spintronics, it offers a unique opportunity 
to pursue novel functionality since spin and momentum 
are tightly related there. In fact, several proposals have 
been made such as the quantized magneto-electric ef- 
fcct^^S, giant spin rotatioiiii, magnetic properties of the 
surface stated, magneto-transport phenomen£feSri2^ ther- 
moelectric properties''^, and superconducting proximity 
effect including Majorana fermionsi^""— . 

Among properties unique to Dirac fermions. Berry 
phase effect is noticeable. Recently, it has been clari- 
fied that Berry phase plays a pivotal role on anomalous 
transport in ferromagnets.^° Now, Berry phase effect is 
recognized as one of fundamental effects and is responsi- 
ble for a wide variety of phenomena.— In topological in- 
sulator. Berry phase effect can manifest itself in physical 
quantities by attaching ferromagnet or making thin-film 
topological insulator—"—. Since magnetization of ferro- 
magnet can be controlled by an applied magnetic field, 
one may expect a remarkable Berry phase effect in mag- 
netic transport on the surface of the topological insula- 
tor. Moreover, topological insulators discovered so far are 
also known as good thermoelectric materials. Therefore, 
the surface of topological insulator offers a good arena to 
study spin caloritronics effect. 

In this paper, we investigate transverse magnetic heat 
transport on the surface of the topological insulator on 
which a ferromagnet is attached. We present general ex- 
pressions of the Peltier and the thermal Hall conductiv- 
ities, which are reduced to simple forms at low temper- 
atures. It is clarified that the Peltier and the thermal 
Hall conductivities show non-monotonous dependence on 
temperature. At low temperature, they have a linear 
dependence on temperature. From the behavior of the 
Peltier conductivity at low temperature, one can esti- 
mate the magnitude of the gap induced by time-reversal 
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FIG. 1: (Color online) Schematic of the model (upper) and 
the dispersion of massive Dirac fermion (lower). 



symmetry breaking. In addition, we find that the Peltier 
conductivity can be used to map the Berry phase struc- 
ture. These results open up possibility to control trans- 
verse thermal currents magnetically. 

We consider transverse heat transport on the surface 
of the topological insulator on which the ferromagnet is 
attached as shown in Fig. [T] Note that time-reversal 
symmetry breaking by the ferromagnet is essential to 
produce Hall effects. Heat current under weak electric 
field E and a small thermal gradient VT is determined 
by 
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Here, we focus on transversal heat transport, namely 
Nernst-Ettinghausen and Leduc-Righi effects, which are 
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characterized by the Peltier and the thermal conductiv- 
ity tensors, a and k, respectively. These quantities can 

be calculated a a^^i^° Uxy = -^f^ci, = —-^^02 where 




Here, fir, Cfcr, ^^'^ f ^-re the Berry curvature with the 
band index r, the dispersion of the surface Dirac fermion, 
the chemical potential, and the Fermi distribution func- 
tion {f{e) = {e^ + respectively. These expressions 
are reduced to: 



irJjSr^^^ [/3(£fc.-M)/(£..-/^)+log(l + e-^(---^))] , (3) 

(4) 



klT f dk 



y + (3^{ekr - A*)'/(£fcr - ^i) - [log(l + g-^^^-^^))] ' - 2Li2(l - fiSkr - A*)) 



where Li2(z) is the polylogarithm function. 
In the limit of T — >■ 0, we have 



and 
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Here, ct^j, is the Hall conductivity. This confirms the 
Wiedemann- Franz law for the Hall currents. We also 
find the relation 



6 OKxy 



(7) 



The coefficient for the T-linear dependence of axy is 
the sum of the Berry curvature at the Fermi energy. 
Therefore, the Peltier conductivity can be used to map 
the Berry phase structure by varying fi. On the other 
hand, the thermal Hall conductivity n^y is proportional 
to the sum of the Berry curvature over the occupied 
states, as is similar to the Hall conductivity (Jxy Note 
that the expressions obtained above are general, applica- 
ble to any nearly degenerate band structures, e.g. narrow 
gap or magnetic semiconductors where the Berry phase 
plays a significant role.~ 

Now, let us consider the following Hamiltonian as a 
general model of Dirac fermion on the surface of the topo- 
logical insulator under the exchange field: 



^ = E4 [(a-kK 



(b • k)ay + (c • k)a,] Cfc (8) 



where a, b and c are real 3D vectors, k = {kx,ky, A)^ , 
and Ux,y,z are Pauli matrices in spin space. Without loss 
of generality, we can set Cz — 1- Now, let us assume two 
fold rotational symmetry ai kx = ky — which imposes 
az — bz = Cx = Cy = 0. The disp ersion is given by 
£kT = r^(a • k)2 + (b ■ k)^ + (c • k)^ with r = ±1 which 
is illustrated in Fig. [TJ After some algebra, we obtain 
the Berry curvature: 



ilr = 



(ax b)^A 



(9) 



This represents the Berry curvature stemming from Dirac 
fermion under the exchange field along 2;-axis. Notice 
that in the presence of the in-plane exchange field, by a 
proper translational operation, we can set Gz = bz = 
in general. With a proper transformation, we can derive 
expressions of the Peltier and the thermal Hall conduc- 
tivities by substituting into Eqs.Q and (|4]) the Berry 
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where k 



\k\ and k represents 2D wavevector: k = 
(kx^ky)* . As is expected, the magnitude of the Berry 
curvature takes its maximum at fc = 0. The extremal 
values of VLt have opposite signs for r = +1 and r = — 1. 



At low temperatures, we obtain a simpler form of a 
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kse [ dk tA , tt kse A 

axy = ^ — T-T jT—T^y -5{h-tV A^ + k^) =T — — -T (11 

" 3 h J (27r)2 2[A2 + A:2]2 12 h ^2 

when the chemical potential is within the band, and axy = otherwise. The thermal Hall conductivity at low 
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FIG. 2: (Color) Peltier conductivity axy in unit of -j^ (left panels), and normalized thermal Hall conductivity K^jy/A in unit 

of (right panels). Thin black lines in the upper panels represent the chemical potential as a function of T for fixed electron 
densities. In lower panels, we show temperature dependence at /^/A = 2 (lower left) and that at /i/A = 0.5 (lower right). The 
straight lines in the lower panels represent the low-temperature asymptotic behavior in Ea. (lll|l and Ea.(|13p. 



temperature is given by 

3 h J (27r)2 ^ 2[A2 + fc2]2 



When the chemical potential is inside the gap, the ther- 
mal Hall conductivity reads 



12 h • 



(13) 



In the left panels of Fig. [21 we show the Peltier con- 
ductivity axy numerically calculated with Eq.([3]) in unit 
of ^j^. The Peltier conductivity axy is an odd function 
of /X. The sign change across — reflects the oppo- 
site sign of the Berry curvature for the upper and lower 
bands. We find a region with large magnitude of axy 
This can be understood as follows. As we noted above, 
the Peltier conductivity axy at low T is proportional to 
the sum of the Berry curvature at the Fermi energy. At 



finite temperatures, thermal broadening changes the S 
function in Ea.(|lip into a distribution over the region 
M ~ £fcr| < aT {a = 0(1))- On the other hand, the Berry 
curvature takes its maximum at = 0. Therefore, when 
—A < fi — aT < A < ^, + aT is satisfied, the contribu- 
tion of Berry curvature from the upper band gives large 
axy On the other hand, if /i — aT < —A < fj, + aT < A 
is satisfied, the contribution of Berry curvature from the 
lower band dominates, and axy becomes a large negative 
value. This is qualitatively consistent with the upper left 
panel of Fig. [51 Thin black lines in the upper panel 
represent chemical potential as a function of T for fixed 
electron densities. Along the lines of fixed electron den- 
sities, nonmonotonic dependence on temperature is seen. 
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Temperature dependence at /x/A = 2 is shown in the 
lower left panel. The straight line in the lower panel 
represents the low-temperature asymptotic behavior in 
Ea.(|lip. Low temperature behavior is well described by 
the analytical expression Ea.([TT|). Thus, from the prefac- 
tor of the temperature, one can estimate the magnitude 
of the gap A. 

The right panels in Figure [5] show normalized thermal 

Hall conductivity Kxy/lS. in unit of numerically calcu- 
lated from Eq. The thermal Hall conductivity n^y is 
an even function of /i. The Kxy becomes large for large /i 
and A. Due to thermal broadening effect, the step func- 
tion would be modified in the region |/i — ekT\ < a-T in 
Eq.(12). As we noted earlier, the thermal Hall conductiv- 
ity Kxy is proportional to the sum of the Berry curvature 
over the occupied states. Since the Berry curvatures for 
the upper and lower bands have opposite signs, they be- 
come small when /x is well below —A or well above A. For 
-A < fi+aT < A or -A < ^-aT < A, the contribution 
of Berry curvature from the lower band dominates over 
that from the upper band. Thus, we expect large thermal 
Hall conductivity in this region. On the other hand, the 
prefactor T in Eq.(12) suppresses the magnitude of K^y 
at low temperature. In this way, we can understand the 
behavior of Kxy in the upper right panel of Fig. [21 From 
the plots, we also find a ~ 3. Thin black lines in the 
upper panel represent chemical potential as a function of 
T for fixed electron densities. Along the lines of fixed 
electron densities, we find non monotonic dependence on 
temperature. Temperature dependence at fi/A. = 0.5 is 
shown in the lower right panel. The straight line in the 



lower panel represents the low-temperature asymptotic 
behavior in Ea.(|13p. which shows good agreement with 
the numerical calculation at low temperature. 

As for experimental realizability, the proximity- 
induced exchange field is expected to be around 5- 
50meVi^ If we set A =10 meV, then the parameter 
range used in this paper is fi = —50 ~ 50 meV and 
T = 1 ^ lOOK, which are within reach of the present-day 
technique. When thermoelectric effects are measured, 
there may be additional contributions to the thermoelec- 
tric effects such as phonon drag. However, these contri- 
butions are negligible at low temperatures. 

In summary, we have investigated transverse magnetic 
heat transport on the surface of topological insulator on 
which a ferromagnet is attached. We have presented gen- 
eral expressions of the Peltier and the thermal Hall con- 
ductivities, which are reduced to simple forms at low tem- 
peratures. It is shown that the Peltier and the thermal 
Hall conductivities show non-monotonous dependence on 
temperature. At low temperature, they have a linear 
dependence on temperature. From the behavior of the 
Peltier conductivity at low temperature, one can esti- 
mate the magnitude of the gap induced by time-reversal 
symmetry breaking. Moreover, we have found that the 
Peltier conductivity represents the Berry curvature at the 
Fermi energy and that it can be used to map the Berry 
phase structure. These results open up possibility to con- 
trol transverse thermal currents magnetically. 
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